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Abstract 

We discuss techniques for obtaining a numerically stable evaluation of the 
fully differential cross section for the virtual photon correction to single hard 
photon bremsstrahlung in two different computational schemes. We also com- 
pare the role of finite mass corrections in these schemes. 

1 Introduction 

Electron-positron annihilation to fermions, e + e~ — ► //, plays a critical role in 
extracting precision electroweak physics from e + e~ colliders. Precise calcula- 
tions are needed for the final data analysis of LEP and the anticipated ILC. 
Thus, the basic process must be augmented by radiative effects, in particular 
bremsstrahlung from a fermion line, e + e~ — > ffj- More photons may be added 
as needed, including virtual photons exchanged in internal lines to the desired 
order. Single hard bremstrahlung is also of importance in "radiative return" 
experiments, where the radiated photon is used to reduce the effective energy 
of the collision, permitting a range of energies to be studied at a fixed-energy 
accelerator. This technique is useful at B factories and high-energy colliders. 

In this note, we will be concerned with the accurate calculation of one-loop 
corrections to single real photon emissions, where there is one extra virtual 
photon in the process. We will discuss some aspects of calculating these 
corrections in a manner stable enough to permit high-precision Monte Carlo 
comparisons to related results. [21 IS] We will also discuss the role of finite- mass 
corrections, and compare our approach to that of Ref. [3], examining the role 
of mass corrections in the level of agreement found for these results. 
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2 Virtual Corrections to Bremsstrahlung 



The one-loop virtual correction to initial or final-state bremsstrahlung was 
calculated in Ref. using helicity spinor methods, which provide an effi- 
cient representation of massless fermion scattering, including a "magic" choice 
for the photon polarization vectors which eliminates many terms from the 
calculation. [31 03 E] The amplitudes were simplified using the symbolic manipu- 
lation FORM,|7] and the scalar one-loop Feynman integrals were evaluated us- 
ing the FF package. |S] Eventually, these integrals were replaced by the analytic 
expressions in Ref. as shown in the appendix of that paper. The amplitudes 
are then evaluated by the /C/C Monte Carlo program, [0] which squares and sums 
them numerically when creating events to obtain the scattering cross section. 

The initial state virtual photon contribution to the cross-section may be 
expressed as 
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with s = (pi + P2Y1 

[pi-Pi) 2 -, U2 = (P2-P3) 2 , ^ = 2pi-k/s, where pi,p 2 ,P3,PA, k are the momenta 
of the e~,e + ,f,f, and 7. The explicit mass corrections in Eq. |21are obtained 
using the method of Ref. as discussed in Sec. 3. The matrix element for 
hard photon initial-state emission with one virtual photon may be expressed as 
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where fi are scalar form factors and 1{ are spinor factors defined in Ref. pp. 

The expressions for the fi include differences between logarithms and dilog- 
arithms with arguments which are very similar in collinear limits, and these 
differences are frequently divided by the collinear factors n, so that the result 
is finite in the collinear limits. Evaluating such expressions in a numerically 
stable manner requires expansions where appropriate. A set of functions use- 
ful for this purpose are the logarithmic and dilogarithmic difference functions 
Li n (x) and Si n (x,y) introduced in Refs. and defined recursively via 



Lf (x) = ln(l + x), 
Lf n+ iO) = - (Lf n (x) - Lf„(0)) , 
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Sfo(z,y) = Sp(x + y), (5) 
Sf n+ i(x, y) = - (Sf n (x, y) - Sf„(x, 0)) , 

with Sp (x) the dilogarithm (Spence function). A set of identities and expan- 
sions for these functions may be found in the appendix of Ref. Thus, 
for example, the form factor /o (the only term in Eq. 0] which survives in the 
collinear limits), may be expressed as, without mass corrections, 

/ \2 \ 2 

/ = 21X1—^ + 3- 27T (L-l-nr)-L 2 -l + — 
V m e J 3 



(1 — ri)(l — r 2 ) I \1 — vj j l-r 2 
(3w + -^M Lfx (-«) + - - R 1 (n,r 2 ) + r 2j Ri(r 2 ,ri) (6) 

with L = ln(s/mg) the "big logarithm" of a leading log expansion, v = r± + r 2 
the fraction of the beam energy radiated into the hard photon, and 



Ri(x,y) = Lfi(-x) jln (^-^) - 2^ 



1 — X \1 — x 1 — X ) 

The parameter A is a photon mass cutoff for the infrared divergence. The 
expression for /q appearing here and in Ref. is analytically identical to the 
expression in Ref. PP, but is preferable for numerical evaluation, because it can 
be evaluated in a stable manner in collinear limits. 



3 Finite Mass Corrections 

It may not be immediately clear that there is any need to consider the finite 
mass of the electron in high energy scattering: for the ILC, m\j s is of order 
1(T 12 . But in fact, in any process where collinear photon emission is possible, 
the electron mass cannot be neglected, regardless of the scattering energy. This 
is because integrating terms of the form m^/(p ■ k) 2 over photon momentum 
k always gives contributions of order 1. Obtaining precise results for collinear 
emission therefore requires mass corrections to be added. Nevertheless, most 
finite mass terms are negligible, so it makes sense to use a massless helicity- 
amplitude technique, but supplement it with a procedure^] to restore the 
essential collinear mass terms. 

For ISR, the net result is to add a mass term to the squared amplitude of 
the form 

M (m) 



E 



(8) 
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where the sum is over the two incoming fermion lines. The effect on the form 
factors is that the spin-averaged value of /o receives an additional mass term 

/f v(m) 2m l {l-v){rl + rl) 

srir 2 (1 - r\Y + (1 ~ r 2) 

x [(/o) + In (J^j [ln(l - v ) - 1] - | ln(l - V ) + ^ ln 2 (l - v) + l| . 

This representation of the mass corrections is very compact and well suited to 
Monte Carlo integration. 



4 Comparison of Virtual Photon Correc- 
tions to ISR 

Another expression for the virtual photon corrections to initial state radiation 
(ISR) has been calculated in Refs. [2J El using a "leptonic tensor" and an ex- 
pansion in factors of m 2 e /p, L ■ k. This result incorporates the same processes 
as Ref. ^ and claims the same order of exactness, but was obtained by very 
different means, so it provides a valuable cross-check. 

An analytic comparison |1 11 I13| showed that the massless parts (without 
explicit electron mass corrections) of the results were identical in collinear limits 
(small pi ■ k), which means that the results are identical to NLL order (order L) 
in an expansion in a leading-log expansion. In this limit, the functions f± and 
/ 2 in Eq. vanish and /o simplifies greatly. The massless spin-averaged NLL 
limit is 

Re(/or" = L^illt-D-l'-ll^ ,^-"' 
\ m£ J 3 1 + (1 — r\Y 

+ Z r2 \ + 2 Inn ln(l - r 2 ) + 21nr 2 ln(l - n) (10) 

1 + (1 - r 2 ) z 

- ln 2 (l - n) - ln 2 (l - r 2 ) + 31n(l - n) + 31n(l - r 2 ), 

and mass corrections can be added via the collinear limits of Eq. 

The mass corrections are more difficult to compare analytically, both be- 
cause the mass corrections of Ref. are much longer, and because those ex- 
pressions include terms proportional to ra\{j>i ■ fc)~ 3 which are absent in Ref. PQ. 
Rewriting the mass corrections using the functions Lf n and Sf n shows that all 
such factors actually cancel, and the two expressions for the mass corrections 
agree to NLL order. [TT] in the sense that all terms producing at least one factor 
of L upon integration agree. In the case of mass corrections, this statement is 
not as strong as saying that the collinear limits agree. In fact, we have shown 
that the soft collinear limits agree analytically, but the L° hard collinear limits 
are not identical. 

Detailed numerical comparisons |llllTllll5j have been made using the /C/CMC, 
both at 1 GeV and 500 GeV, corresponding to typical energies for B factories 
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Figure 1: Differences between the exact and massless NLL results for the virtual correction to 
single photon 1SR, both with and without mass corrections, in a /C/CMC run of 10 8 events. Results 
are shown for two different energy scales, yfs = 1 GeV and 500 GeV. JMWY is the expression of 
Ref. and KR is the expression of Ref. |3| , and ti max is the limit on the radiated photon energy, 
relative to the beam energy. 



and the ILC, respectively. Since the massless NLL limits are known to agree, 
we have subtracted the NLL limit (Eq. IKJj) in each case, and plotted the re- 
maining contribution, both with and without explicit mass terms. The results 
of MC runs for 10 8 events are shown in Fig. 1, normalized to the Born cross 
section. 

The plots show that the massless results agree to within 10~ 6 for all but 
the last bin of the 500 GeV run. The main difference is evident in the mass 
corrections. The maximum difference between the integrated cross sections is 
found to be 2.2 x 10~ 5 at 1 GeV and 1.6 x 10 -5 at 500 GeV, in units of the 
Born cross section. 1 While this is an excellent level of agreement, it is clear 
that there is some disagreement in the mass corrections, which are responsible 
for almost the entire difference. It would be desirable to understand the origin 
of this difference more completely. 



Acknowledgments 

This work was supported in part by U.S. Department of Energy grant DE- 
FG02-05ER41399. 



For radiative return, it is more relevant to compare to the hard ISR cross-section oi, which is 
0.113 cr B orn at 1 GeV and 0.980 cr Bo m at 500 GeV. 
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